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The accuracy of equation-of-state formulations is compared for theoretical total energies or 
experimental pressure- volume measurements for H2, Ne, Pt, and Ta. This spans the entire range of 
compression found for minerals and volatiles in the Earth. The Vinet equation is found to be most 
accurate. The origin of the behavior of different equation-of-state formulations is discussed. It is 
shown that subtle phase transitions can be detected by examining the residuals from an equation-of- 
state fit. A change in the electronic structure of Ta is found at high pressures using this procedure, 
and a possible new transition in H2. 



I. INTRODUCTION 



In geophysics and high-pressure research, experimental or theoretical data consisting of pressure, temperature, 
and volume triplets PVT are parameterized to a functional form for ease of interpolation and extrapolation. These 
equations of state are then used to compute phase diagrams, or are used in geodynamic or shock compression modeling, 
etc. A recent book by Anderson (1995) comprehensively reviews equation-of-state formulations, and comprehensive 
reviews and comparisons of equations of state are given by Hama and Suito (1996), Stacey et al. (1981), and Duffy 
and Wang (1998). The primary purpose of this paper is to discuss the reasons for the accuracy of the most-used 
formulation, based on the Birch equation of state (Birch 1978) and compare with the Vinet equation of state (Vinet 
et al. 1987). Jeanloz (1988) has previously shown that the Birch equation of state and the Vinet equation of state can 
be similar up to moderate compressions. Here we show that this breaks down at high compressions, and for highly 
compressible materials the Vinet equation of state is considerably more accurate. A logarithmic equation of state 
(Poirer and Tarantola 1998) and another exponential equation of state (Holzapfel 1996) are discussed and compared 
with the Vinet and Birch equations. Finally, we show how subtle transitions can be detected by deviations from an 
equation of state. 

The above equations of state are appropriate for isothermal data. To include thermal expansivity there are three 
main approaches that have been used. The first is to fit isotherms, and then tabulate or fit the parameters Vq, Kq, 
and Kq as functions of temperature. Since experimental data are not always collected along isotherms, that method 
is most amenable to analysis of theoretical results. Secondly, one can assume Debye and Mie-Gruneisen theory, and 
fit parameters for O^), 7, and q to include thermal effects. Thirdly, one can directly model the thermal pressure. See 
Anderson (1995), Hama and Suito (1996), and Duffy and Wang (1998) for further discussion of thermal equations of 
state. 

In order to test equation-of-state formulations, it is important to study large compressions, because most common 
equations of state will work reasonably over small compression ranges. Thus we choose to study the very compressible 
materials, hydrogen and neon. We also consider tantalum and platinum, which are useful pressure standards. 



II. ISOTHERMAL OR STATIC EQUATIONS OF STATE 

The Birch equation of state (Birch 1978) is based on a series expansion of pressure 



P(/) = 3ifo/(l + 2/)^ 



in terms of the Eulerian strain /, where 



1+2(^0-4)/ + 



(1) 



/ 



H- 



(2) 



and Vo, Kq, and Kq are the zero pressure volume, bulk modulus, and bulk modulus pressure derivative. Truncated 
as written it is called the "third-order" Birch equation. The fourth-order Birch equation, which includes another term 
in and is useful for theoretical equations of state, but when applied to experimental data, the added parameter 
usually induces severe correlations in parameters so that they lose physical significance. 
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The Vinet equation (Vinet et al. 1987) is derived from a scaled approximate form for the energy: 



E{r) = -A£;(l + a*)exp| 



r - To 



(3) 



where AE is the binding energy, and r is the length per electron. This gives: 



P{x) = 3Ko {I - x) x-^ eiip 

1 

where x — (^^^ . The energy can then be expressed as 
4KoVo 



- [K', - 1) (1 - x) 



(4) 



E = Eo 



2VoXo(/<o - 1)" (5 + iK{x - 1) - 3a;) exp(--(i^^ - \){x - 1)). 



The Holzapfel equation of state (Holzapfel 1996) is similarly given by 

P{x) = iKax-^il - x) exp [{cx + Co) (1 - x)] 



(5) 



(6) 



where cq and c are chosen to give K' and the limiting Fermi gas behavior as a; — > 0. If c = 0, one gets a 3 parameter 
(Vq, Kq, and Kq) equation of state that behaves better at extreme compression (Hama and Suito 1996): 



P{x) = SKqx '^{1 — x) exp [co {1 — X 

"3 



(7) 



3Kqx '^(1 ~ x) exp 



{K^~3){l-x) 



Finally we consider the logarithmic equation of state recently proposed by Poirer and Tarantola (1998), derived 
nilarly to 
third-order. 



similarly to the Birch equation, but where the strain is defined as e = In^^. rather than as in cq. giving at 



P = Ko 



lny + (^^)(ln-) 



(8) 



Now we examine the application of these equations of state for several materials. First we consider the extremely 
compressible behavior of hydrogen. Hemley et al. (1990) showed in analyses of x-ray diffraction data for hydrogen 
that the Vinet equation of state is considerably more accurate than the Birch equation. The experimental equation 
of state, (room temperature data is corrected to 4 K), is shown in Figure 1 along with equation-of-state fits, and the 
fitted parameters and are shown in Table ||. For most of the fits, Vb was fixed at its known value, 23.0 cm'^/mole 
(Silvera 1980). The best fit is given by the Vinet equation of state, with acceptable fits using Birch or Holzapfel, but 
the logarithmic equation of state fails completely. When the parameter Vq is relaxed for the Vinet equation there is 
no significant change in the fit, and the parameters become less well determined. Figure lb shows the residuals of 
the fits (the logarithmic residuals go off scale as shown). Interestingly, there is a peak in the residuals at 40 GPa, 
which may be a subtle transition, or change in compression mechanism, previously undiscovered, or a small problem 
with a subset of the data. Note that this deviation is completely invisible on the P-V curve itself. This illustrates 
the usefulness of having a good equation-of-state formulation; even small deviations may indicate possible transitions 
that should be examined more closely. Secondly, it illustrates that one only expects equations of state to work well 
under compression with no underlying transitions. If there are even subtle electronic or structural transitions, these 
will affect the equation of state significantly. 

One of the tests of an equation-of-state fit is the accuracy with which it reproduces known zero-pressure parame- 
ters. To this end, high-pressure static compression data can be compared with low-pressure elasticity measurements, 
including data from ultrasonic experiments, Brillouin scattering, or in some cases high-precision static compression 
data. Useful low-pressure static compression data exist for hydrogen because of its very high compressibility. Swenson 
and Anderson (1974) reported Kqt = 0.17±0.06 GPa for n-H2 from volumetric strain measurements to 2.5 GPa at 4.2 
K. The results are in excellent agreement with those of Wanner and Meyer (1973), who obtained Kq = 0.174±0.010 
GPa for single-crystal n-H2 at 4.2 K using ultrasonic methods. The latter number provides the best comparison 
with the present analysis (e.g., for n-Il2). These results are also close to zero-pressure Brillouin scattering results 
obtained at T=4 K for p-Hs by Thomas et al. (1978); the latter obtained Kqs = 0.173±0.001, with Kqt = 0.162 GPa 
(corrected for isothermal conditions). Udovidchenko and Manzhclli (1970) performed accurate static compression 
(volumetric) measurements on P-H2 from to 18 MPa (down to 6 K) and obtained Kqt — 0.186±0.006 GPa. The 
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neutron diffraction study of P-H2 to 2.5 GPa by Ishmaev et al. (Ishmaev et al. 1983) gave Kqt = 0.186±0.03 and 
K'^rp = 6.33±0.2. The Vinet Kq is close to that determined directly, but the Birch and Holzapfel Kq's are way too 
high. The deviation in i^o.T with the Vinet equation may be due to being thrown off by the equation of state glitch at 
40 GPa, or may be due to insufficient flexibility in the Vinet equation over this large compression range. In order to 
examine this possibility, we considered the extended Vinet equation given by Moriarty (1995) and Vinet et al. (1989): 



P = 3Ka (1 - x) exp 



-(K^-l){l-x)+P{l-xf + j{l-xf 



(9) 



where P = (-19 + 18/^0 + 9Kl^^ + 36KqKI{) and fit the H2 data varying /3 as well as Vq, Kq, and Kg. As shown 
in table |, Kq = 0.15 ± 0.02, now in excellent agreement with the ultrasonic bulk modulus. The peak in residuals 
survives this higher-order fit. Other extensions to the Vinet equation are possible, and it would seem more physical to 
add further terms to the prefactor of eq. ^, such as 6(a*)~'^, rather than terms in the exponential, but such extended 
equations of state will not be considered further here. 

Secondly, we consider another soft material, Ne, from fitting theoretical Linearized Augmented Plane Wave 
(LAPW) (Singh 1994) local density approximation (LDA)(Hedin and Lundqvist 1971) results. Figure 2 shows the 
fitted equation of state, the residuals of the fit, and the resulting PV curves, equation-of-state parameters are shown 
in Table |l[ Only the Vinet and Birch equations are compared since the logarithmic equation of state fails above for 
H2, and the energy expression obtained from integrating the Holzapfel PV equation is not closed form. The Vinet 
equation is superior to the Birch equation again, and does an excellent job simultaneously over the large energy scale 
of large compressions, and the small energy scale near the minimum. We also show results for the fourth-order Birch 
equation, which in spite of an additional parameter, is still not as good as the Vinet equation. As more parameters 
are added, the quality of the fit improves, but the correlations among indi vidu al parameters increases. 



Next we consider equations of state of two stiff metals, Ta and Pt (tables |lVj and III and Figures 3 and 4). Energies 
were computed using the LAPW method and the Generalized gradient approximation [GGA-PBE, (Perdew et al. 
1996)]. The Vinet equation fits better than the Birch equation for Pt, but for Ta the Birch equation fits better. 
However, examination of the residuals for Ta (Fig. 4c) shows that there is a peak in the residuals at ~10.5 A'^, and 
closer study shows a change in the occupied bands (and thus the Fermi surface) around that volume, where the t2g 
states at F move below the Fermi level, Ep, (Fig. 5) so the improved fit using Birch is a side effect of this anomaly. 
Fitting only the data above 13 A'^ shows a greatly improved fit, and a better fit for Vinet than Birch. This again 
illustrates the utility of examining residuals from well founded equations of state such as the Vinet equation in order 
to find subtle transitions. 



III. DISCUSSION 



We have considered a range of materials from very compressible to more incompressible, covering the entire range of 
compression of minerals and volatiles in the Earth. An equation of state reflects the underlying interactions potential 
among the ions and electrons that make up a crystal. Thus, a useful way to compare different equations of state is 
to compare the implied assumptions about the interaction potential. First we consider the Birch equation in some 
detail. As normally expressed, it is not obvious that the Birch equation in terms of Eulerian strain 



E = 



N ''^ ri 



n=0 



(10) 



— 2ti 

where 77 = is actually a series for energy E in (^^^ ^ . For example, if = 3 



111/113 

E = ao--ai + -02 - gfls + I -^ai - -02 + ^03 ) 77 



1 



1 



(11) 



The I powers arise from distance squared in finite strain theory, so one can think of the Birch equation as a series 
in "inverse length squared." The reason only even terms in length are included is in order to preserve rotational 
invariance in the strain energy expansion for general strains. Thus, the Birch equation assumes that the underlying 
potential can be represented as a series in (1/r^"). The commonly used third-order Birch equation (with parameters 
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Vq, Kq, and Kq) includes n = 1,2, and 3. Now, it is well known that polynomials can behave poorly (wiggle) and 
extrapolate poorly (go crazy outside the range of the fit) , and this is in fact a problem with using the Birch equation 
to too high an order. Note also that it is not a convergent series (unless the coefficients obey ^ < /i which is not 

expected in a fit). Under compression, rj — V/Vq < 1 so that for larger n the terms r]~^ blow up. The truncated series 
may still do well, as it does up to moderate compressions, if the truncated series represents the interatomic potential 
accurately. However, there is no fundamental reason to expect (l/r^") to well represent interatomic interactions, so 
it is not surprising that the Birch equation is not perfect. Indeed, the energy expression from the Birch equation has 
a non-physical hump at expanded volumes that then decays with increasing volume, whereas the Vinet and Holzapfel 
equations converge smoothly to a constant at large volumes, being consistent with physically based potentials. The 
logarithmic equation of state unphysically diverges with expansion. 

Another question that often arises with regard to the Birch equation is why Kq — 4 often works so well. Note that 
for moderate compressions (say 10% or less) higher order terms in eq. |l| are quite small. For example, if 77 = 0.9, then 
/ = 0.036, so that deviations from Kq — 4 are only refiected in a few percent or less of the pressure, and higher order 
terms, with factor = 0.0013 are only affected at the 0.1% level. 

The Lagrangian strain 



1 f i 



1 



(12) 



can be used instead of the Eulerian strain (Thomsen 1970), but the Lagrangian equation of state is not satisfactory 
at high pressures since it saturates with increasing density. The Lagrangian strain gives a series in E 



E 



N 

E 

n=0 



N 



n=0 



- 1 



N 



(13) 



which has the advantage that the series converges in 77^ for 77 < 1, but since ^ r^" does not look like an interatomic 
potential, it does not work well compared with the Eulerian strain. 
The Hencky (logarithmic) strain (Poirer and Tarantola 1998) 



dl 

den = —;eH 



ml 



(14) 



shares a similar shortcoming. Though it looks reasonable as a measure of finite strain the effective potential, ^ (Inr)" 
is not very physical. The Hencky potential can be expanded to look like — a(l/r) + fc(l/r)('^^''~^^ (Poirer and Tarantola 
1998). The attractive part is long-ranged and ill-conditioned and the repulsivepart is too soft. 

In contrast, the Vinet equation is based on a potential at the outset (eq. y). This potential was introduced by 
Rydberg (1932), who used it not for an equation of state (Stacey et al. 1981), but as a form for the miramolecular 
potential in H2 and other simple molecules for obtaining solutions to the anharmonic Schrodinger equation for fitting 
and comparison to molecular vibrational spectra. The Vinet equation of state works surprisingly well for a wide range 
of types of materials, and for compressions of up to 77 = 0.1, as has been discussed in detail (Hama and Suito 1996; 
Vinet et al. 1987). Note the sign errors in eq. 102 of Stacey et al. (1981). 

One failure of the Vinet equation is that it does not give the proper behavior at even greater compressions, and 
does not merge into electron gas (Thomas-Fermi), or quantum-statistical (Kalitkin and Kuz'mina 1972) behavior at 
extreme compression. Holzapfel found expressions that do have the proper limiting form, such as eqs. |^ and |^ above. 
The underlying energy expression, however, of even the simplest Holzapfel equation is exceedingly complex 



E{x) = lx-^Koi9exp 
"3 



(K'q - 3) (1 - x) 



+3 exp 



{K'o-3)x 



{-4 + 2{-5 + 3K[,) X 
3 



(15) 



(15 - UK^ + 3K'^) x^ Ei 



[K'^-3)x 



where 



The fact that the Holzapfel equation works less well than Vinet for H2 and other materials (Hama and Suito 1996) 
up to almost 10-fold compression suggests that this underlying potential docs not accurately represent the interaction 
potential in solids, though it has the correct limiting behavior at extreme compressions. 

An equation of state that agrees with the quantum statistical model (Kalitkin and Kuz'mina 1972) at extreme 
compression, which is obtained by assuming homogeneous compression of spherical atoms and expanding the Hartree- 
Fock equations in density and dropping terms of order and higher, was developed by Hama and Suito (1996). This 
should become accurate at lower pressures than Thomas-Fermi. We did not test this here, but in Hama and Suito's 
tests, it seems to work well from to 10 fold compression or more fore rare gases, simple metals, and ionic compounds. 
Further exploration of this approach is warranted for extreme compressions. Their equation of state is: 



where x = (y-j , as above, and ^ is determined from the quantum statistical model (Kalitkin and Kuz'mina 1972). 



Finally, in comparing equations of state, note that correlations among Vq, Kq and Kq are much lower for Vinet 
than Birch, so parameters are better determined. In the Vinet equation, Kq is inside exponential, and Kq in the 
prefactor, whereas in the Birch equation Kq and K'q appear as products. In any case, it is very difficult to resolve Vq 
from V-P data, especially for a high-pressure phase unstable at P = 0. Data usually only extends to P = at best, 
and it is thus difficult to obtain good values of Kq and Kq (if not for a robust equation-of-state formalism) since one 
is trying to find first and second derivatives from one-sided data. This is one advantage of theory; results are available 
at negative pressure and Vq can be well constrained. 

For strains lc;ss than 30%, it probably doesn't matter what equation of state you use, as was emphasized by Jcanloz 
(1988), but parameters will still be better determined with the Vinet equation (see also Hemley et al. 1990). For large 
strains, the Vinet equation is best, and forms such as the Holzapfel equation are required at extreme compressions. 
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TABLE I. Fitted equation-of-state parameters for H2 pressure-volume experimental data. 



X 



Vq (cm^/mole) 



Kq (GPa 







Vinet 
Vinot 
Birch 3 
Birch 4 
Birch 4 
Holzapfel 
Logarithmic 
Ext. Vinet 



0.1901 

0.1895 

0.6554 

0.2356 

0.4033 

0.4126 

13.98 

0.1623 



23.0 

22.4±2.6 

23.0 

23.0 

23.0 

23.0 

23.0 

23.0 



0.250±0.004 

0.279±0.15 

0.694±0.009 

0.28±0.02 

0.162 

0.47±0.01 

0.08±0.36 

0.15±0.02 



6.56±0.02 

6.49±0.35 

3.990±0.007 

5.7±0.2 

7.77±0.03 

4.95±0.03 

64.3±292.4 

8.0±0.4 



-9.4±1.2 
-23.4±0.3 



-21. ±2. 



TABLE II. Fitted equation-of-state parameters for Ne LAPW (LDA) energy (Ryd, 1 Ryd=13.605 eV))-volume results. 



Vinet 



Birch 3 



Birch 4 



I. 35x10-' 
8.598±0.074 

II. 2±0.8 
6.07±0.14 

Mn.2)" 

-2o(j,7(i."KS I (1.00(12 



Vo (cm"^/molc) 
Ko (GPa) 
K 

KoA',' 

Eo (liyil) 



1.02x10-" 
8.641±0.008 
8.94±0.06 
7.192±0.009 

M6.0)" 

-2o0.7oii:! I (1,(10( 



3.03x10-" 
8.592±0.013 
9.24±0.18 
7.51±0.11 

-21.1+0,98 
-2o0.7(i:.l I 0.00003 



*Not varied, result of formulation. 



TABLE III. Fitted equation-of-state parameters for Pt LAPW (GGA-PBE) energy (Ryd)-volume results (V>10A^). 



Vinet 



Birch 



exp. 



Vo (cm^/mole) 
Ko (GPa) 

Eo (Ryd) 



1.47x10-' 
9.448±0.007 
248.9±0.7 
5.43±0.02 

-36893.5648±0.0002 



3.51x10-' 

9.445±0.011 

238.6±1.1 

5.474±0.03 

-36893.5641±0.0002 



9.09 
278 
5.6 



'■(Holmes et al. 1989) 



TABLE IV. Fitted equation-of-state parameters for Ta LAPW (GGA-PBE) energy (Ryd)-volume results. 



Vinci 



Birch 



Vinet (V>13A'') 



Birch (V>13A^) 



exp. 



Vq (cm' / mole) 
Ko (GPa) 
Ko' 

Eo (Ryd) 



9.32x10-' 
11.026±0.009 
188. ±1.8 
4.08±0.04 

-31252.3337± 0.0004 



2.71x10-' 
11.058±0.005 
188. ±0.9 
3.81±0.01 

-31252.3339±0.0002 



8.08x10-" 
11.057±0.002 
192.3±0.2 
3.82±0.02 

-31252. 33425±0. 00004 



I. 31x10"" 

II. 062±0.03 10.865 
190. ±0.3 195±5 
3.77±0.02 3.4±0.1 
-31252.3342±0.00005 



'■(Cynn and Yoo 1999) 
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IV. FIGURE CAPTIONS 



Figure 1. Equation-of-state of hydrogen, (a) Experimental P vs. V. and equation-of-state fits, (b) Residuals of 

equation-of-statc fits. Note that the logarithmic equation-of-state residuals go off scale. 

Figure 2. Equation of state of neon from LAPW total energy computations within the LDA. (a) Energy versus 
volume and third-order Birch and Vinet fits. The inset (inset area is shaded in the larger figure) shows the that the 
Birch equation does not match the data at low pressures, whereas the Vinet works at both low and high compressions, 
(b) Pressure volume relations from the fits shown in (a) as well as a fourth-order Birch fit. In spite of the fact that 
aU fits look good to the eye at high pressures in figure (a), the derived pressures differ by over 150 GPa at extreme 
compression, (c) Residuals for the fits. Note the large and systematic deviations in the third-order Birch fit. The 
fourth-order Birch is almost as good as Vinet, but with one additional fitting parameter. 

Figure 3. Static equation of state for Pt to extreme compressions computed using LAPW and GGA. (a) Energy 
versus volume and Vinet, third- and fourth-order Birch fits. The "low P" Vinet fit is to volumes greater than 10 
A^. (b) Pressure versus volume from the fits. The Vinet and fourth-order Birch give similar results over the extreme 
pressiire range. A Vinet fit to a smaller pressure range gives better constrained equation of state parameters and 
smaller residuals, but the pressures are very close to those computed using the large data range, (c) Residuals of the 
Vinet fits, and the lower pressure third-order Birch fit. Note that there are no systematic deviations. 

Figure 4. Static equation of state for Ta computed using LAPW and GGA. (a) Energy versus volume. Note that 
to the eye both fits appear perfect, (b) Pressure versus volume. In spite of the apparently excellent fits in (a), the 
pressure is significantly different under extreme compression, (c) Residuals of the fits. An electronic transition was 
hidden under the apparently smooth equation of state. There is a peak in the residuals at 10.5 A^. When data is fit 
at V>13A^, the fit greatly improves and the residuals are much smaller. 

Figure 5. Band structures of Ta at (a) V=12.66 A^ (5 GPa) and (b) V=9.3 A^ (460 GPa). There is a major change 
in the occupied states as the t2g states move below the Fermi level Ej? at high pressures, crossing at a volume of 11.67 
A^ (200 GPa). This change is responsible for a glitch in the equation of state. 
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